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Abstract 

We describe an integrable model, related to the Gaudin magnet, and its re- 
lation to the matrix model of Brezin, Itzykson, Parisi and Zuber. Relation is 
based on Bethe ansatz for the integrable model and its interpretation using 
orthogonal polynomials and saddle point approximation. Large N limit of 
the matrix model corresponds to the thermodynamic limit of the integrable 
system. In this limit (functional) Bethe ansatz is the same as the generating 
function for correlators of the matrix models. 
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1. Introduction 

Matrix models in the large A^-limit have been studied from different points 
of view since the seminal papers of t'Hooft [1] and Brezin, Itzykson, Parisi 
and Zuber [2] . There is no doubt of their importance in 2D gravity and string 
theory (see e.g. [3], [4], [5] and [6] for reviews on their different aspects). 
Importance of CFT and of integrable systems for matrix models have also 
been known for some time (see e.g. [4], [5]). Most recently there is a great 
interest in integrable spin chains and Bethe ansatz in relation with large N- 
hmit of = 4 SYM (see e.g. [7], [8], [9] and references therein). This possibly 
may help to understand more precisely the large N conjecture of Maldacena 
[10]. 

The purpose of this paper is to make precise the relation between the large 
A/'-limit of the one-matrix model and Bethe ansatz for some integrable model 
related to the classical rational r-matrix (Gaudin type model). This could 
provide a toy model of the above mentioned relation between the Heisenberg 
model and A/" = 4 SYM. 

We start with a very short review of the one-matrix model following [3]. 
Two basic methods of its solution, the saddle point approximation and the 
method of orthogonal polynomials, contain implicitly all information about 
some particular integrable model. Then we recall the version of Bethe ansatz 
appropriate for models related to classical r-matriccs [11], [12], [13]. This 
can be viewed as some specific limit [11] of the quantum /2- matrix Bethe 
ansatz [14], [11], [15], [16] which is appropriate for Heisenberg chain type of 
models. Finally we will go on with an explicit description of the relevant 
integrable model and discuss algebraic as well as functional Bethe ansatze 
for it. We will relate in a precise way the large iV- limit of the matrix model 
and the thermodynamic limit of our integrable system. In this limit Bethe 
eigenfunction give the generating function for correlators of the matrix model. 
Also the norm of the Bethe eigenstate is nicely related to the large N expansion 
the partition function Zjy. 



2. Matrix model of BIPZ 

One-matrix model is treated in detail in many places [3], [4]. So we can be 
very brief. The task is to perform the integral over N x N hermitian matrix 
M 

Zn^ j d""" Mexp [-{N/g)TrViM)] , (1) 

where V{M) is a general polynomial potential V{M) = + J2k>3(^kX'^ and 
the integration is over independent real variables which are the real and 
imaginary parts of matrix elements of M. The standard strategy is to diago- 
nalize the matrix M, i.e. write M = ?7diag(Ai, Aat)?/^ with an unitary U. 
Noticing that the integrand depends only on N eigenvalues A^ and integrating 
out the unitary group one gets 

N 

Zn= I l[d\exp [-{N/g)^V{X,) + J^HK - A,)] . (2) 

1=1 i i^j 
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The logarithmic part in the exponent justifies the saddle point approximation 
for N ^ oo. This gives the configuration with the minimal energy 

|Ex4^-nA.). (3) 

In order to solve equation (3), it is convenient to introduce the resolvent 
function 

i ^ ' 

One also introduces the averaged density of eigenvalues p(A) 

i 

If we assume for simplicity that the potential V has only one well (or at least 
one deepest well), then the support of p(A) will be in some interval [a, 6] (in 
general it can be a collection of disjoint intervals). We have relations 

W(z) = fdX^^ (6) 



and 



The solution for W{z) is 



m^)^-L/,.mv(-»)(-^), (8) 

^ ' 4m J z-x ^{x-a){x-b) ^ ' 

where the integration contour encloses the cut [a,b]. Finally the solution to 
(3) is given by the density (7). See the above cited references and [17] for 
more details concerning the solution of this equation. 

An alternative method of solving (1) is to use orthogonal polynomials 
-Pn(A) = A" + ... with respect to the measure 

exp[-{N/g)V{X)]dX. (9) 

Partition function is then simply expressed in terms of their norms squared 

Zjv = 7V!so...s^-i. (10) 

This follows easily from the integral representation for P„ which will be given 
later (37), (38). See again [3], [4], [17] for a more explicit method of deter- 
mination of these polynomials. This is based on recurrence relations between 
orthogonal polynomials which can be solved in the large A^-limit. This lit- 
tle knowledge about the one-matrix model is enough to understand it as an 
quantum integrable model when N is large. 



3. Bethe ansatz 

Before describing our integrable model we give a brief review of the relevant 
Bethe ansatz following [13]. Let L(A)bea2x2 matrix 

[C{X) -A{X)) 

its matrix elements being operators depending on complex parameter A, acting 
in some Hilbert space Ti. We assume the following form of the commutation 
relations for the matrix elements of L(\) 

[L{X) ® /, J ® + [r(A - 12), L(A) ® / + J ® = , (12) 

where / is the unit 2x2 matrix and r(A — fi) is the classical r-matrix 

/l 0\ 
^0010 
A 1 

\0 ly 

Note that the commutator in (12) is a matrix commutator respecting the 
operator nature of matrix elements of L{X). Further, we put 

T(A) = iTrL(A)^ (14) 

where Tr is the "2 x 2 matrix trace", as a consequence of (12) we have the 
following useful commutation relations 

[A(A), S(/.)] = -Y^B(^i) + Y^B{X) , (15) 

[BiX),B{f2)] = 0, (16) 
[B{X),C{X)] = 2A'{X) (17) 



(13) 



and 



[T(A),S(/.)] = Ai^,)B{X)-- B{^,)A{X) 

A — fJ, A — ylt 

^ .A{X)B{i,) + -^B{X)A{i,). (18) 



X — ji X — fj, 

Operator T(A) is the generating function of integrals of motion in involution 
as 

[T(A),r(/.)] = 0. (19) 
Next we assume that there is the so-called pseudo- vacuum |0) e 7i such that 

C(A)|0) = 0, (20) 
A(A)|0) = a(A)|0) (21) 

and hence 

r(A)|0) = (a(A)2 + a'(A))|0) (22) 
holds true. Algebraic Bethe ansatz (ABA) is then given by a family of vectors 

\Xi,...,Xn)^B{Xi)...B(Xn)\0). (23) 
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These are, as easy verified, with help of (12-18), eigenstates of T{X) with 
eigenvalues t(A) 

r(X,^-2^^^^^^^^^ia(Xf + aW (24) 
if and only if the numbers {Ai, Ajv} satisfy Bethe conditions 

^^,^a{\)^Y.j:^-^- (25) 



The norm squared is equal to 

difi 



(Ai, AatIAi, Xn) — det 



(26) 



If we introduce the polynomial g(A) = (A — Ai)...(A — Ajv) equation (25) can 
be rewritten as differential equation 

g" - 2ag' + (a^ - a')q = rq . (27) 

4. The model and its relation to BIPZ at large N 

Comparing Bethe conditions (25) and the minimal energy condition (3) we 
see the first hint that there might be an integrable model hidden in the large 
A^-limit of the matrix model. To describe this integrable system we need to 
give a representation of the algebra of (12). Let us fix an arbitrary integer 
K and consider the space of functions of K variables a;_x}. Our 

operators ^(A), -B(A) and C{X) will be polynomials of orders K, K — 1 and 
K — 1, respectively, acting on this space 

-K-l -K -K 

A{X) = ^nA-"-^ , B{X) = SnX-""-' , C(A) = Y CnA-"-^ , 



n=— 1 n=— 1 



with coefficients 



(28) 



Bn^Xn, (30) 

Cn — x^dm—n—idi + 2a^9^_„ . (31) 

Here the summation rule is assumed and sums run over all integers from — oo 
to K, assuming 9„ = and = whenever n > —1 and {a_i, ..a_x-i} 
are arbitrary real constants. Note that A_k-i = a_K-i- we will assume this 
to be nonzero. What we have here is actually the free field representation of 

a factor algebra of sl{2) with c — 0. ABA is apphed directly to our model. 
Obviously pseudo- vacuum is the constant function |0) = 1 and we have 

-K-l 

a{X) = Yl ■ (32) 

n=-l 

To our integrable model also Sklyanin's functional Bethe ansatz (FBA) is 
applicable. FBA in its full generality doesn't need a pseudo-vacuum. In our 
case it gives nicely the separation of variables. We will not repeat here the 



general discussion of [12] and just state the result. Since B{X) is a polynomial 
of order K — 1, 

5(A) = x_i^(A - Mi)...(A - , (33) 

x^K and its zeros {ui, ...,uk-i} are independent as functions of variables 
{x-i, X-k}- We perform the change of variables 

X-k} ^ {ui, ...,uk-i, X-k} ■ (34) 

If q{u) = {u — Xi)...{u — Ajv) is a polynomial solution of order N to (27) then 
the product 

i){ui, ...,uk-i) = x^K q{ui)...q{uK--i) (35) 

is an eigenfunction of T(A) and the corresponding eigenvalue is given by (24). 
The relation between the two version of Bethe ansatz is pretty obvious in our 
case. We go form (23) to (35) as follows 

N K-1 

I Ai, = s(Ai)...s(A^).i = n n - ^^■) 

i=i j=i 

K-l N K-l 

j=i i=i j=i 

^(-lfi^~')x'[j,^l^(u,,...,UK-i). (36) 

I both cases (algebraic or functional) the Bethe condition (25) just makes 
sure that poles of t(A) in points {Ai,...,Aiv} vanish. This must happen as 
the operator T(A) doesn't have poles in these points. If we put a_i = 0, 
a_2 = N/g and a_„ = for n > 3, so that we have a = j^V, then Bethe 

conditions (25) and the minimum energy condition (3) are identical. 

We can go even further. Following Szego's book on orthogonal polynomials 
(and being inspired by matrix model) we consider a system of orthonormal 
polynomials Pn{u) with respect to the measure (9) given as 

(£>._i£>J-V^ 

Pn{u) = j X (37) 

77/. 

/n 
Y\ dtti {u - ai)...{u - an) exp [- ^{N/g)V{ai) + ^ ln(Q;j - aj)] , 
1=1 i i^j 

where 

^ „ n+l 

= u^TTv. / n ^^^^ [- J2W9)y{(^i) + H'^i - «^-)]- (38) 

^ ' i=l i i^j 

Obviously we do have 



Pn{u) 



Dr 



.1/2 



D. 



n-l 



Pn{u) (39) 



and (10) follows trivially. Let us now assume N ^ oo. Again the logarithmic 
term in the potential allows to apply the saddle point method to the integrals 
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defining Pn{u), which are generating functions of correlators for matrix modeL 
This is done straightforwardly with the result 

Pn{u) = {u- X,)...{u - A^) + 0{[N/g])-' . (40) 

Here {Ai, .... X^} must again minimize the energy, hence they satisfy (3). Re- 
calling the functional Bethe ansatz (35) we see that the polynomials 

Pn{ui)...Pn{uk-i) q{ui)...q{uK^i) (41) 

in the large iV-limit and apart from an inessential factor x^j^ they give Bethe 
eigenvectors ip{ui, ...,uk-i) of T(A). Also (see (26)) 

N/2 

(Ai...Aiv|Ai...AAr)-^/' X 

exp [-{N/g) V{Xi) + H\ " A,)] (1 + 0{g/N)). (42) 

i i^j 

Although Bethe states (23) or (35) are eigenstates of our integrable model 
for any N here we recovered them only in the large A'^-limit. For our rep- 
resentation labels {a_2, a_A--i} this means that they grow proportionally 
to N as this goes to infinity. At the same time N is the number of vacuum 
excitations. Hence we can interpret the large A^-hmit as the thermodynamic 
limit of our integrable system. Clearly in this limit the polynomial q{u) is the 
generating function of correlators for the matrix model. 
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